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ABSTRACT 

We report the preliminary results of a survey for water vapor in a sample of eight C stars with large 
mid-IR continuum fluxes: V384 Per, CIT 6, V Hya, Y CVn, IRAS 15194-5115, V Cyg, S Cep, and 
IRC+40540. This survey, performed using the HIFI instrument on board the Herschel Space Obser- 
vatory, entailed observations of the lowest transitions of both ortho- and para-water: the 556.936 GHz 
lio — loi and 1113.343 GHz In — Ooo transitions, respectively. Water vapor was unequivocally de- 
tected in all eight of the target stars. Prior to this survey, IRC-l-10216 was the only carbon-rich AGB 
star from which thermal water emissions had been discovered, in that case with the use of the Sub- 
millimeter Wave Astronomy Satellite (SWAS). Our results indicate that IRC-l-10216 is not unusual, 
except insofar as its proximity to Earth leads to a large line flux that was detectable with SWAS. The 
water spectral line widths are typically similar to those of CO rotational lines, arguing against the 
vaporization of a Kuiper belt analog (Ford & Neufeld 2001) being the general explanation for water 
vapor in carbon-rich ACB stars. There is no apparent correlation between the ratio of the integrated 
water line fluxes to the 6.3 micron continuum flux - a ratio which measures the water outflow rate - 
and the total mass-loss rate for the stars in our sample. 

Subject headings: circumstcllar matter — stars: AGB and post-AGB — stars: abundances — subniil- 
limeter: stars 



1. INTRODUCTION 

Over the past two decades, observations of the en- 
velopes of AGB (asymptotic giant branch) stars have re- 
vealed a molecular composition that is still not entirely 
understood. In early models for the chemistry within 
the circumstellar envelopes of AGB stars, the inner en- 
velope was assumed to possess a composition determined 
by thermochemical equilibrium (TE) within the stellar 
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photosphere, while the outer envelope, exposed to the in- 
terstellar radiation field, exhibits a time-dependent gas- 
phase chemistry that is driven by the photodissociative 
effects of ultraviolet radiation (e.g. Glassgold 1996). In 
TE, the carbon-to-oxygen ratio is crucial in determin- 
ing the photospheric composition; cool oxygen-rich stars, 
with C/0 < 1, are expected to have photospheres that 
are dominated by CO and H2O, while those of carbon- 
rich stars will be dominated by CO and C2H2. However, 
a variety of observations have revealed anomalous abun- 
dances for several molecules that cannot be explained by 
models in which the inner envelope entirely reflects the 
photospheric abundances in TE and the ensuing chem- 
istry proceeds in the gas phase. More recent models in- 
volving shocks (Willacy & Cherchneff 1998; Cherchneff 
2006), grain surface reactions (e.g. Willacy 2004), the va- 
porization of orbiting objects (Ford & Neufeld 2001), or 
the penetration of UV radiation through a clumpy out- 
flow to the inner envelope (Decin et al. 2010) have been 
variously invoked to explain the observed abundances of 
such species as CO, HON, SiO, CS, H2CO, OH and H2O 
in the envelopes of both 0-rich and C-rich AGB stars. 

While most of the principal circumstellar molecules 
are easily detected by means of ground-based observa- 
tions, water has been more elusive. Although water 
mascr emissions from oxygen-rich stars can be detected 
from the ground, thermal water emissions are detectable 
only with satellite observatories. As expected, the In- 
frared Space Observatory (ISO) detected luminous water 
emissions from a variety of oxygen-rich AGB stars (e.g. 
Barlow 1999, Neufeld et al. 1999). More surprisingly, 
later observations performed (Melnick et al. 2001) with 
the Submillimeter Wave Astronomy Satellite (SWAS) 
- and subsequently confirmed with the ODIN satellite 
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TABLE 1 

Sample of C stars surveyed for water vapor 



Star 




R.A. 




Dec. 


D 


M 


^sys 


Date (2010) 


Duration (s) 


Noise (: 


mK) f 




(J2000) 


(J2000) 


(pc) 


(a) 


(b) 


557 


1113 


557 


1113 


557 


1113 


















GHz 


GHz 


GHz 


GHz 


GHz 


GHz 


V384 Per 


03 


26 29.5 


+47 


31 50 


740 


5.5 


-16.2 


Sep 1 


Sep 2 


3972"^ 


987 


2.4 


14 


IRC+10216'= 


09 


47 57.4 


+13 


16 44 


110 


15 


-25.5 


May 4,11 


May 11 


3150"^ 


3200" 


2.7 


14 


CIT 6" 


10 


16 02.2 


+30 


34 19 


320 


1.9 


-1.6 


May 11 


May 11 


4719 


987 


1.9 


16 


V Hya 


10 


51 37.2 


-21 


15 DO 


600 


2.1 


-17.0 


Jun 23 


Jun 9 


3972** 


987 


3.4 


15 


Y CVn 


12 


45 07.8 


+45 


26 25 


220 


0.15 


21.1 


Jun 23 


Jun 9 


3972"^ 


987 


2.2 


16 


IRAS 15194-5115 


15 


23 04.9 


-51 


25 59 


590 


18 


-15.0 


Aug 2 


N/A<= 


3972** 


N/A-^ 


2.1 


N/A" 


V Cyg 


20 


41 18.3 


+48 


08 29 


530 


2.5 


14.3 


Jun 23 


Mar 5 


3972'* 


987 


2.6 


16 


S Cep 


21 


35 12.8 


+78 


37 28 


720 


3.4 


-15.5 




May 12 


N/A<= 


987 




13 


IRC+40540'^ 


23 


34 27.7 


+43 


33 02 


780 


7.0 


-17.0 


Jun 23 


Jun 9 


3972<* 


987 


2.5 


15 



^ Mass loss rate in units of 10 ^ Mq yr ^ 

^ Systemic velocity in km s~"^ relative to the local standard of rest 

IRC+10216, CIT6, and IRC+40540 are known also as CW Leo, RW LMi, and LP And respectively. 
^ Divided into two observations of equal duration with slightly different LO settings (see text) 
^ N/A — not available; planned observation not yet carried out 
^ R.m.s. in a 1 km s^^ bandwidth 



(Hasegawa et al. 2006) - have led to the discovery of wa.- 
ter vapor in the envelope of the nearest carbon-rich AGB 
star with a massive envelope, IRC-l-10216. Prior to the 
recent launch of the Herschel Space Observatory (Pil- 
bratt et al. 2010), this detection of the lowest, lio — loi; 
rotational transition of ortho-water at 556.936 GHz by 
SWAS and ODIN remained the only case in which ther- 
mal water vapor emission had been detected from a C 
star. 

The Herschel Space Observatory presents an opportu- 
nity to search for water vapor in C stars with unprece- 
dented sensitivity; with a primary mirror of collecting 
area ~ 30 times that of the SWAS mirror, and with its 
Heterodyne Instrument for the Far Infrared (HIFI; de 
Graauw et al. 2010) employing cryogenic mixers of much 
lower noise temperature than those flown on SWAS, Her- 
schel makes it feasible to search for water in AGB stars 
considerably more distant than IRC+10216. In addition, 
Herschel provides access to multiple transitions of water 
vapor. As part of the HIFISTARS Key Program, we car- 
ried out a survey for water vapor in a sample comprising 
eight C stars in addition to IRC+10216 - V384 Per, CIT 
6, V Hya, Y CVn, IRAS 15194-5115, V Cyg, S Cep, 
and IRC+40540 - with the goal of determining whether 
IRC+10216 is unusual (in addition to being nearby), 
or whether water vapor is widely detectable in the en- 
velopes of C stars. Our survey involves observations of 
the lowest transitions of both ortho- and para-water, the 
556.936 GHz lio - loi and 1113.343 GHz In - Oqo tran- 
sitions respectively, toward a selection of C stars with 
large mid-IR continuum fluxes. Table 1 lists the sources 
in our sample, together with estimates of their distances, 
D, mass-loss rates, M, and systemic velocities relative to 
the local standard of rest. Although the quantity M/D^ 
is fairly well constrained, uncertainties in D make both 
parameters quite uncertain - a fact that is demonstrated 
by the wide dispersion of estimates given in the literature 
- but in the interest of uniformity we adopt the values 
tabulated by Groencwegcn et al. (2002), supplemented 
by those of Schoicr & Olafsson (2001) for the case of Y 
CVn (for which estimates were not given by Groenewc- 
gen et al. 2002). With the exception of Y CVn, which is 



classified as a J-type variable carbon star of spectral type 
CV5J (with V indicating variability and J indicating a 
large ^■^C/^^C ratio), all the stars in the sample are of 
spectral type CV6 or CV7. V Hya is also of a somewhat 
different character than the other sources; it is a semi- 
regular variable (SRa), quite probably a binary (Sahai 
et al. 2003), and possesses a high velocity outflow and 
bipolar symmetry that suggests it is in an early stage of 
post-AGB evolution (Kahane et al. 1996, Knapp et al. 
1997). The mass- loss rate given in Table 1 is for the 
component with outflow speed ^ 15kms^^ (identified as 
the "normal slow wind" by Knapp et al. 1997). 

In this Letter, we report preliminary results of this sur- 
vey, based on the observations performed to date toward 
the eight target stars. The detection of the In— Oqo water 
transition toward V Cygni in the very first observation 
performed in this survey has been reported previously 
(Neufeld et al. 2010). The data acquisition and reduc- 
tion are described in §2, and the results are presented 
in §3. In §4, the results are discussed with reference to 
various scenarios that have been proposed for the origin 
for water vapor around C stars. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations reported here were carried out during 
the period March - September 2010 as part of the HIFIS- 
TARS Key Program. We used the HIFI instrument in 
dual beam switch (DBS) mode to observe two rotational 
transitions of water vapor towards each source listed in 
Table 1. The 556.936 GHz lio - loi and 1113.343 GHz 
111 — Ooo transitions were targeted respectively in lower 
and upper sidebands of Bands lb and 4b of the HIFI 
instrument. The telescope beam was centered at the co- 
ordinates listed in Table 1, with the reference positions 
located at offsets of 3' on either side of each source. The 
beam sizes were 38" and 19" (half power beam width) 
respectively for the 556.936 GHz and 1113.343 GHz tran- 
sitions. The data were acquired using the Wide Band 
Spectrometer (WBS), which provides a spectral reso- 
lution of 1.1 MHz and a simultaneous bandwidth of 
~ 4 GHz. For each observation, the date and total du- 
ration are given in Table 1, together with the r.m.s. 
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Fig. 1. — Black histogram: Hcrschel/HIFI spectra of the lio — loi 
556.936 GHz transition of ortho-H20 (left) and the In — Ooo 
1113.343 GHz transition of para-H20 (right) observed toward V384 
Per, IRC-l-10216, CIT6, and V Hya. Blue histogram: analogous 
spectra for CO J = 10 — 9, where available, scaled to the same 
peak antenna temperature. The spectra have been continuum- 
subtracted and rebinned to a spectral resolution of 1 km s~^. 
Doppler velocities are expressed in km s~^ relative to the systemic 
velocity of the source, and antenna temperatures (vertical axis) are 
shown in mK. Dotted lines indicate expansion velocities from the 
literature ( Groenewegen et al. 2002). 
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Fig. 2.— Same as Fig. 1, for Y CVn, IRAS-15194, V Cyg, S Cep 
and IRC-l-40540. For Y CVn and IRAS-15194 the blue histogram 
shows CO J = 6 — 5 instead of J = 10 — 9 



noise achieved in a 1 km bandwidth. As discussed 
in Neufeld et al. (2011), the data were processed using 
the standard HIFI pipeline to Level 2, providing fully 
calibrated spectra with the intensities expressed as an- 
tenna temperature and the frequencies in the frame of 
the Local Standard of Rest (LSR). For several observa- 
tions, we used a pair of local oscillator (LO) frequencies, 
separated by a small offset, to confirm the assignment 
of any observed spectral feature to either the upper or 
lower sideband of the (double side band) HIFI receivers. 
The resultant spectra were coadded so as to recover the 
signal-to-noise ratio that would have been obtained at a 
single LO setting. Spectra obtained for the horizontal 
and vertical polarizations were found to be very similar 
in their appearance and noise characteristics and were 
likewise coadded. 

3. RESULTS 

Figures 1 and 2 show the continuum-subtracted spec- 
tra obtained from each of the observations listed in Ta- 
ble 1, with the water spectra plotted in black. To im- 
prove the signal-to-noise ratio, each spectrum has been 
rebinned to a spectral resolution of 1 km s~^. In each 
case, the velocity scale is relative to the systemic veloc- 
ity of the source (Table 1, as given by Knapp et al. 1998 
for Y CVn and Groenwegen et al. 2002 for the other 
stars.) The blue histograms show, for comparison, the 
CO J = 6 - 5 ( in the case of IRAS 15194-5115 and Y 
CVn) or CO J = 10 — 9 (other stars) spectra obtained 
with HIFI, scaled to have the same peak antenna tem- 
perature as the water spectra. No HIFI spectra of CO 
are yet available for V384 Per. Water vapor is clearly 
detected in each of the eight target stars by means of 
one or both of the observed transitions. As predicted 
by excitation models that we performed in anticipation 
of Herschel (Gonzalez- Alfonso, Neufeld, & Melnick 2007; 
hereafter GNM) , and given the sensitivity of Herschel as 
a function of frequency, the lio — loi transition provides 
the most sensitive probe of water in C stars (although in 
the single case of S Cep, the detection of water vapor was 
obtained through observations of the In — Oqo transition 
alone, the lio — loi having not yet been observed.) In 
Table 2, we list the integrated antenna temperatures and 
line fluxes for each transition, together with the 6.3//m 
continuum flux obtained from the ISO observations of 
Sloan et al. (2003) or - in the cases of CIT6, V Hya and 
IRAS 15194-5115 - estimated from the S.Sfim fluxes re- 
ported by Monnier, Geballc & Danchi (1998). The abso- 
lute flux calibration for ISO spectroscopy was accurate to 
within ^ 30%. The continuum flux is of particular inter- 
est because theoretical models for the excitation of water 
in C stars (Agiindcz & Cernicharo 2006; GNM) indicate 
that the observed lio — loi and In — Ooo rotational emis- 
sions are excited primarily by radiative pumping through 
the 1^2 vibrational band at 6.3 /iin. 

4. DISCUSSION 

Although there are real differences between the H2O 
and CO line profiles for some of the spectra shown in 
Figures 1 and 2 (most notably for the IRC-|~10216 spec- 
tra which have the highest signal-to- noise ratio) , the line 
widths are typically quite similar. In particular, the half 
width at zero intensity is typically close to literature es- 
timates of the expansion velocity (outer dotted vertical 
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TABLE 2 
Water line intensities and fluxes 



otar 


/ TAdv 


Line flux" 


^6.3 


■- 

X 




( K km 


s-1 ) 


(lO-i'^ W m-2) 


(lU Jy) 


15 




557 


1113 


557 


1113 




C 

iz 




GHz 


GHz 


GHz 


GHz 




15 
Z5 
C 


V384 Per 


0.64 


1.12 


0.55 


1.95 


0.46 




IRC+10216 


10.3 


14.6 


8.70 


25.22 


22.3 


o 

U 


CIT 6 


2.17 


3.81 


1.86 


6.62 


4.00 


C 


V Hya 


0.28 


0.65 


0.24 


1.13 


1.70 


O 
U 


Y CVn 


0.11 


0.27 


0.09 


0.47 


0.37 


'e 


IRAS 15194-5115 


1.12 


N/A'' 


0.96 


N/A'' 


1.00 




V Cyg 


1.10 


1.74 


0.94 


3.02 


1.04 


CD 


S Cep 


N/A'' 


2.54 


N/A'' 


4.41 


0.94 


X 
15 


IRC+40540 


0.29 


0.37 


0.25 


0.64 


0.40 


OJ 


^ We conservatively estimate the typical flux 


uncertainty 


as better 




than 15%, except for 


V Hya, 


Y CVn, 


and lRC+40540. 


In these 


"O 
OJ 


sources, the 1113 GHz 


flux uncertainties arc dominated by 


noise; we 


O 


estimate them as 20% 


30%, an 


d 25% (1 


cr) res 


pcctivcly. 




cn 

03 


^ Planned observation 


not yet 


carried out 






c 



lines in each panel). This feature of the spectra indi- 
cates that water-emitting gas is typically present over a 
wide range of outflow directions, including streamlines 
directed close to the linc-of-sight, and argues strongly 
against the vaporization of a Kuiper belt analog (Ford 
& Neufeld 2001) being the general explanation for water 
vapor in carbon-rich AGB stars. Any collection of orbit- 
ing icy objects with a flattened structure (like the Kuiper 
belt) would release water that would be entrained in the 
outflow and carried along streamlines with small incli- 
nations to the equatorial plane of the structure. Viewed 
nearly edge-on, such an anisotropic distribution of wa- 
ter vapor could indeed yield a spectral line of width 
comparable to that of the CO lines; but viewed nearly 
face-on, the water lines would be considerably narrower. 
Thus, while the vaporization of a Kuiper belt might ac- 
count for the water line profile observed in any individual 
source, the statistical ensemble now available from our 
survey is not consistent with Kuiper belt vaporization as 
a general explanation for water vapor. (Note, however, 
that based upon multitransition observations available 
for IRC-l-10216 but not for the other sources in our sam- 
ple, Decin et al. 2010, and subsequently Neufeld et al. 
2011, have argued that the distribution of water vapor 
in that source is inconsistent with the vaporization of 
small icy objects on circular orbits being the source of 
the water in that source.) 

Detailed calculations by GNM predict luminosities for 
the lio — loi a-iid In — Oqo transitions that arc roughly 
proportional to the 6.3 y^tm continuum luminosity, and 
that the ratio of integrated line flux to 6.3/im contin- 
uum flux is an increasing function of the water mass-loss 
rate. Because coUisional excitation is negligible except in 
the densest CSEs, the line fluxes show almost no depen- 
dence upon the gas density; thus the ratio of integrated 
line flux, F;, to continuum flux, Fg.a, is nearly indepen- 
dent of the total mass- loss rate. In Figure 3, we plot 
the Fj/Fg.a ratio for each star as a function of the total 
mass-loss rate, M. Here, we used the results tabulated 
in Table 2, with integrated line fluxes in units of 10^^^ 
W m^^, and the 6.3 /xm continuum fluxes are in units 
of 10^ Jy; thus the plotted ratio has units of MHz. The 
plotted values of M are the literature estimates tabu- 
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Fig. 3. — Ratio of integrated line flux to 6.3/im continuum flux for 
the 557 GHz (black diamonds) and 1113 GHz (red asterisks) water 
lines, as a function of the estimated total mass-loss rate (Table 1) 

lated in Table 1, with error bars indicating an estimated 
uncertainty of 0.2 dex. If the water vapor present in 
these CSEs is produced by chemical means within the 
outflow with some universal H2O/II2 ratio, the Fj/Fg.s 
ratio (measuring water outflow rate) and the total mass- 
loss rate would show a positive correlation. There is no 
discernable correlation in Figure 3, suggesting that stars 
with larger mass outflow rates tend to have smaller wa- 
ter abundances (although it remains unclear whether this 
tendency indicates a causative relationship). Quantita- 
tive estimates of the water outflow rates and abundances 
will require detailed modeling that will be the subject of 
a future paper (and that we hope to constrain further 
with future Herschel observations.) For the two sources 
that have been studied in detail to date, V Cyg (Neufeld 
et al. 2010) and IRC-hl0216 (GNM; Neufeld et al. 2011), 
the derived water outflow rates support the lack of posi- 
tive correlation suggested by Figure 3. While IRC-M0216 
has a total mass-loss rate ~ 10 times that of V Cyg, it 
shows a water outflow rate that is only ^ one-third as 
large; the inferred II2O/H2 abundance ratio in the CSE 
of V Cyg (2 X 10~^) is therefore ~ 25 times larger than 
that inferred for IRC-M0216 (8 x IQ-^). 

HIFI has been designed and built by a consortium of in- 
stitutes and university departments from across Europe, 
Canada and the United States under the leadership of 
SRON Netherlands Institute for Space Research, Gronin- 
gen, The Netherlands and with major contributions 
from Germany, France and the US. Consortium mem- 
bers are: Canada: CSA, U. Waterloo; France: CESR, 
LAB, LERMA, IRAM; Germany: KOSMA, MPIfR, 
MPS; Ireland, NUI Maynooth; Italy: ASI, IFSTINAF, 
Osservatorio Astrofisico di Arcetri- INAF; Netherlands: 
SRON, TUD; Poland: CAMK, CBK; Spain: Observato- 
rio Astronomico Nacional (IGN) , Centro de Astrobiologa 
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(CSIC-INTA). Sweden: Chalmers University of Tech- 
nology - MC2, RSS & CARD; Onsala Space Observa- 
tory; Swedish National Space Board, Stockholm Univer- 
sity - Stockholm Observatory; Switzerland: ETH Zurich, 
FHNW; USA: Cahech, JPL, NHSC. 
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